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STATE-VARIABLE ANALYSIS OF NON-LINEAR 
CIRCUITS WITH A DESK COMPUTER 

INTRODUCTION 

This work was prompted by the need to analyze the transient performance of the 
power regulator unit for the multi-mission modular spacecraft, in particular under the condi- 
tion of short-circuit failure of the switching transi$tor of a power module. Due to the non- 
linearities introduced by the filter inductors and by the solar array characteristics, the state- 
variable programs available for desk computers could not be used. The purpose of this 
work is to fill up this void and to be general enough to handle most nonlinear circuit or 
system analyses. 

The nonlinearities considered here are not restricted to any particular circuit element. 
They may arise from any passive or active source. Wliat tlie program needs is tlie funda- 
mental relationship governing each nonlinearity in tlie form of points on a curve; for example, 
the flux linkage-current relationship of a nonlinear inductance or tlie voltage-current relation- 

I 

’"Ship of a source. 

The starting point of tlie program is a set of flnt-order differential equations and 
algebraic equations describing the system. That is provided by the user. It was therefore 
deemed useful to include in tills document tlie metliodology of writing equations directly 
from a simple examination of a circuit. Examples have been included, where appropriate, 
in order to illustrate such methodology. 

Hie program is interactive and offen many options to tlie user, among which is plotting 
of the results. It was used very successfully for tlie transient analysis of tlie power module 
mentioned above. 

Circuits As Graphs 

The graph representation of a circuit enables focusing on the manner in which the various 
elements of tlie circuit are interconnected. Each node of tlie circuit has its counterpart in 
rile graoh, a node being rite point at which two or more circuit elements join. For tlie 
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piupQwt of this onilysts, all sourcM (voltaic or current) arc considered by themselves, dis* 
sodated from any other circuit elements; besides, their values are presumed to be variable. 

Any passive circuit element (resistance, inductance or capacitance) found connected in 
parallel with a voltage source or in series witli a current source may be removed as it has 
no bearing on the analysis at lumd. The graph of tlie circuit is then obtained by representing 
each circuit element, except those removed, by a line (edge) joining tlie nodes at the terminals 
of that element. This is illustrated in Fig. I . 



Figure 1 . Circuit and its graph, (a) Original ciraut. (b) Modified circuit where elements 7 
and 8 are dropped because they were in parallel with a voltage source and series with a cuirent 
source, respecti^iy. (c) Circuit graph witli 5 nodes and 6 edges. 


The tree of a graph is obtained, using tlie following procedure. Starting out with 
the node configuration only, enough edges are subsequently added to interconnect the 
nodes witliout forming any cl^d patlis. The edges forming a tree are called tree 
branches and the remaining edges are called links. It is, in general, possible to derive a 
large variety of trees from a graph as illustrated in Fig. 2. In (b) the tree branches are 
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Pi|ure 2. Graph and some of its trees, (a) Original 
graph, (b) A (c) Possible, but not all, trees. 


1, 2, 3, and 4, while the links are S and 6. In (c) 2, 3, 5, and 6 are tree branches, 1 and 
4 are links, 

Due to KirchhofTs voltage law (KVL), tlie tree branch voltages may be considered 
independent variables. The addition of any link to a tree produces a unique closed path. 
This enables the vtiting of a link voltage in terms of a unique combination of tree branch 
voltages. Thus, in die tree of Pig. 2(b), the voltage of link 5 con be expressed in terms 
of the voltages of branches 1, 2, 3, and 4, while the voltage of link 6 depends on brancli 
voltages 3 and 4. In Fig. 2(c), link 1 needs branches 2, S, and 6 for its voltage, and 
link 4 branches 3 and 6. Each link can tlius be associated with a unique set of tree 
branches, namely, those branches which lie along tlie closed path defined by the link, 
and form a “tie set." See Pig. 3(a). 

Likewiw, each tree branch may be associated witii a unique set of links to form a 
“cut set" as follows. A closed surface can be found, crossed by that tree branch alone, 
such that the tree nodes are split into two distinct groups. The links crossing that surface 
are those associated ninth the tree branch. Por tlie tree of Pig, 2(b), the following associa- 
tions hold, as illiutrated in Pig. 3(b). 

Tree branch Cut set Link Tie set 

1 1,5 5 1, 2, 3,4,5 

2 2,5 6 3,4,6 

3 3,5,6 

4 4,5,6 
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For tht tiM of Fig. 2(c), tlie iMocii^'ioni are: 

Tree branch Cut set 

2 1,2 

3 3,4 



f 


Unk Tie set 



Figure 3. Use of a tree in deflning: (a) link voitage, 
(b) tree branch current. 


The Concept of State 

At; any instant of tune, tlie amount of energy stored in an energy-storing element is an 
indication of the “state” of tltat eiement. In a capacitance tlie energy stored is C v^, 

»,e 

in an inductance it is 'A Li^. Tlius, tlie voltage v may be considered the state of a 
capacitance, describing completely its present, independent of its past. Similarly, the 
current i may be considered the str,te of an inductance. 

The transition Of an eiement from one state to anotlier requires the flow of energy 
into or out of tlie element, i.e. a certain behavior of tlie electrical quantity which i$ 
not descriptive of its state. In tlie case of a capacitance, the transition from a voitage v ^ 
to a voitage v j necessitates a flow of current whose behavior in time is responsible 
for taking tlie capacitance from the first state to tlie second. An inductance undergoes 
a transition from a state i ^ , to a state i 2 on tlie heeis of a voltage performance in time. 

From a different standpoint, it can be said that a knowledge of the behavior of 
the current of a capacitance between two instants of time t^ and t 2 is not enough to 
determine the voitage (state) of die capacitance at time t 2 ; it is essential to know also 
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tl!« State from which it started, i .e. tlie voltefc at t| . The same holds true for an inductance 
where the voltage variations across it between times t| and t 2 cannot determine its state 
i] unless its state i i is known. As for a resistance, no particular behavior in time is 
needed of eitlier Of its electrical quantities, i and v . It is a memoryleu element with 
no need for a concept of transition. The idea of state is tlierefore foreign to it. Let 
us move on now to a more complex configuration, the circuit. 

Tlie state of a circuit may be conceived as the set of states of all of its energy-storing 
elements. The electrical quantities, which ore capable of effecting a transition of tiie circuit 
from one state to anotlier (the capacitance currents and tlie inductance voltages) are inter- 
related by the circuit topology and tlie source values. Those electrical quantities may be 
obtained at any time if the state of the circuit and the source values are known. As an 
example, the circuit of Fig. 4(a) is in tlie following state at t ■ 0,1 second: 

Voltage across C • 3 V 

Current through L ■ 0.2 A 



(•) (b) 

Fig. 4. An example: (a) the circuit at any time; 

(b) the circuit at t ■ 0.1 s. 


This state and the source values are shown in Fig. 4(b). Using simple dc circuit analysis 

techniques, the following results are obtained: 

i 1 ■ 0.1 1 A 

i 2 - 0.39 A 

V* ■ 5.0 V 

V b| ■ 4.6 V 
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Hm ctpicitanct currant and the inductance voitafe at t ■ 0.1 1 an thus kndti^, namely, 
0.39 A and 4.6 V, respectively. The transition to the “next" circuit state can then 
presumably be determined. 

In the equivalent circuit of Fig. 4<b), no distinction is made, or is necessary, between 
the d.l V of tKe source and the 3 V of the capacitance state, or between the 0.48 A of 
the source and tlie 0.2 A of the inductance state. It is therefore convenient to include 
tlte source variables, in tlie set deflning tlie circuit state, as the augmented set permits the 
detemUnation of the voltage and current of every element in the circuit. 

State»Variable Approach to Analysis 

In che sequel two assumptions are implicitly made for any circuit; 

(a) no tie set contains only capacitances; (b) no cut set contains only inductances. 

As a result, it is possible to And a tree where all the circuit capacitances are included, 
and none of the circuit inductances. To simplify tiiis first analysis, resistances are assumed 
to be all tree branches or all links. Voltage sources are to belong to the tree, while 
current sources have to be links. The state variables are tlten defined as tlie capacitance 
voltages, the inductance currents, and tlie source variables (voltage or current). 

If the circuit is linear, the next step is to write a set of nnt*order differential 
equations for the state variables, obtained directly from KVL and Kirhoffs current law (KCL). 

For n state variables, denoted by Xj, ..., X„, the rth equation is of the form; 

^ * atl *r2 + - + a»n 

niere are n such equations, which in matrix form may be expressed as: 

X - 

The detailed procedure now follows. First, if the resistances are tree branches, express 
their individual currents in terms of statewariable currents. The voltage is then the current 
expression multiplied by R. On tlie other hand, if tlie resistances are links, their individual 
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voltifM an txpitaatd in tenni of itat*>vari«bl« voitaft*. curront if thin thi voitofi 


•xpiMiion dividid by R. In iithtr enm both tlw volttfe itnd cumnt of tach mistanct 
may be exprantd in tenni of state vaiiabiei only. And now tiie state equations. 

tach capacitance is a tree bnmch and its voltafs v is a state variable. The derivative 
dv/dt is equal to the capacitance cuirent divided by C. Ujiii tlie gnplt method described 
above, tlie capacitance current may be expressed in terms of link currents, I.e, state* 
variable currents only as sought. The same holds true for the inductances. 

Each inductance is a link and its current i is a state variable. Tlie derivative di/dt is 
equal to tlie inductance voltage divided by L. Using the graph metliod, the inductance 
voltage may be expressed in terms of tree*branch voltages, i .s state-variable voltages only. 
Finally, a (set oO first-order differential equation(s) is obtained for each source, obviously 
in terms of its own state variables) only. For more details concerning sources, see 
Hewlett-Packard manual for model 30 calculaton, entitled *‘State Variables PAC,*’ pages 


36-37, 


Example 1. 




The sources are dc with values A and B. A tree is formed including ail capacitances, 
excluding all inductances, and including all resistances. The tree also includes the voltage 
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iouret, but not th« cuirint lourc*. All ititt vaiiablc poUritiw choitn arbitreilly. 


Reaiatance Rj 

cumnt 

I2 


voltafs 

R(i2 

Realatanee R2 

current 

h 


voitafe 


Reiiatance Rj 

currant 

i2 * tj 


voltsfe 

R| (it * ij) 

Resistance R4 

currant 

-ij + i$ 


vottage 

R4 Wj + is) 




I V4 - V| - Rj O2 + ij)J 


<^^3 ■ |R4 Wi + I5) - V| - Rj (ij + ij) -Rji}! 

L2 



£L- ■ 0 

ar^ 


In matrix form, 


• «i 

Vl 


"0 

I/C 

1/C 

0 

0 

vr 

*2 


-l/L, 

KRi+R3)/Li 

-Rj/L, 

l/L, 

0 

i: 

‘3 

m 

-I/L 2 

-R 3 /L 2 

-(R4+R,+R2)/L2 

0 

R 4 /L 2 

'3 

V 4 


0 

0 

0 

0 

G 

V4 

is 

m 0 


0 

0 

0 

0 

0 

is 


coefficient matrix 



Tht initial itita victor ii 

.V ,(0) 

ijCO) whin V |(0) is the initial voitafi of C| 

i](0) ijCO) is thi initial cumnt of L| 

A i](0) is thi initial current of L2 

B 


Example 2. _ 

R| L 



tree 


Here oil resistances are links. 

Resistance Rj voltage ■ 

cumnt ■ (V4 - v,)/R| 

Resistance R; voltage ■ v ( 

cumnt ■ V1/R2 

Resistance Rj voltage ■ V2 

current ■ V2/R3 
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dvt 

dt 

dVi 

dt 

dV4 

dt 

dij 

dt 

dXj 

dt 


J. +‘51 

% 

1 Iv, -v,l 

TT 


-•V4 

-w*U 


for voltage source 


for current source, since 
dis/dt ■ - <u B sin (cut 0) 
d*l5/dt* ■ -w^Bcos (cjt + 


In matrix form, 


dt 


M 

V, 

m 

is 

V4 

‘5 

Xs 


-(l/R, + l/RjVCi 
0 

1 /L 

0 

0 

0 


0 

-1/R3C2 

-I/L 

0 

0 

0 


-1/Cl 

l/c, 

0 

0 

0 

0 


1/C, 

0 

0 

-a 

0 

0 


0 

I/C2 

0 

0 

0 

•cj2 


Tlie initial state vector is 
v,(0) 


where v ,(0) is the initial voltage of C, 

V 2(0) is tlie initial voltage of C2 

13(0) is the initial current of L 


^2 

is 

V4 

‘5 

Xs 


V2(0) 
is(0) 

A 

B cos 0 
-cuB sin 0 

I 

Concerning the status of resistances in a graph, it is not always possible to have them 
all either tree branches or links. A mix is the norm ratlier than the exception. In this 
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caat coniider temporarily tliat the voltages of tree^branch resistances and the currents of 
link resistances ore state variables, Write an equation for the current of each tree*branch 
reaistance and the voltage of each link resistance, using state variables only. Solve for 
the resistance state variables in terms of the remaining state variables, 


Example 3 , 





Let the temporary state variables be Vj for R,, Ig for Rj, v^ for R3, and Vg for R4. 


vs/Ri 

isR: 

^/Rj 


' Vg 4. V, 

' h 

Solve for Vj, i^, V;, and Vg. 


8/R4 


Vj - R,l2 
Rjig - Vg» V, 

v^ ■ Rjij 

Vs/Rs + "«2 -»3 


Vs ■ Riij _ 

, 1 R 4 . R 4 

ig « V, + », i, 

Ra + R4 Ra+ R4 Rj + R4 

V? ■ Mi 

R 

V .1 V. 4- ^2^* U - ^2^4 i. 


Ra+R 4 


Ra ^R 4 


Ra'^R 4 


II 



Thn pnviout proctdutt ia now followad, 


RMutance Ri 


Rtaisti^ice Rj 


Resistance Rj 


Resistance R^ 


dvi «--- 

It C 

5 k 


dt 


Li 

1 


cuirent ■ {2 

voltofe ■ R,ij 

curreiit *» ij 

voltage ■ R]>3 


' iqnc '' ' ■RT+T'7 ’ * ‘‘ " ’’ 

* 

-RiU + *4 - + "■ 




dii ■ — - R 4 „ . R}R4 i R}R4 I , w d i 

^ Lj +R^'» 


dV 4 ■ 0 
dt 


Notice that in all tlie state equations written so for tlie differentiated variables were tlie 
state variables. Tliis is always true of linear systems, not necessarily true of nonlinear 
systems. For the latter a distinction is made between the differentiated variables and 
the state variables. The state equations take the form ^ ■ AX instead of X ■ AX . 

The basic equation characterizing an mductance is v « dX/dt, where X is Uie coil flu,x 
linkage. If X is produced by a coil current i, tlien tlie inductance t ■ X/i. L is not 
constant unless X is directly proportional to i, in which case v ■ dX/dt • d(Li)dt ■ L di/dt. 
For a magnetic core whose B*H curve is available, the X - i relationsliip is derived hi the 
following manner. 



X - N# ■ (NA) B , 
i - (l/N) H 

where N is the number of turns of the coil, 

A is the cross-section of tite mSgnetic core in 

£ is tile mean lengtli of the magnetic core in m. 

’( 

In setting up tlie state equation of an inductance, write dX/dt in terms of tiie state 
variables using KVL. For a nonlinear inductance, tlie differentiated variable is thus X, 
and the state variable i. 

The above expressions for X and i are based on SI units, i.e. B is in tesla (weber/m^ ) 

and H in ampere/meter. Tlie B-H curve is often available in cgs units. In that case, 

X - (NA X 10-«) B 

i - l8/(0.4irN)J H 

where A is in square centimeter, 

8 is in centimeter, 

B is in gauss, 

H is in oersted, 

As for a capacitance, the basic equationriri » dq/dt, where q is the charge in 
coulomb on tlie plates of the capacitance. Tlie charge is related to voltage across tlie 
capacitance by q > Cv, where C is the capacitance in farad. If C is constant (linear 
capacitance), then i ■ C dv/dt; otlierwise, the lelationsliip between q and v must be 
known. In setting up tlie state equation, write dq/dt in terms of the state variables 
using KCL. For a nonlinear capacitance, tlie differentiated variable is thus q and tlie 
state variable v . 

For a resistance v ■ Ri. If R is not constant (nonlinear resistance), one of its 
electrical quantities is taken as a state variable. If the resistance is a tree branch, its 
state variable is the voltage. If the resistance is a link, its state variable is tlie current. 

In either case, tlie otlier electrical quantity is expressed in terms of the state variables 
of the whole problem, using KVL in case of a link, KCL in case of a tree branch. Tliis 
other quantity is obviously not differentiated. 
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Exampit 4, 

Consider Uie circuit of example 1 in ttiis section, Let C, L; , and R4 be nonlinear, 
Using the same tree, assign voltage v^ to R4, V( is a state variable since R4 is a tree 
branch. The state equations become: 

iSl • + Ij 

dt 

““Ri ij '^'^4 (I2 +13) 

dt 

IV 4 - V| - Rj (ij + is) ” R 2 *j 1 

dt L2 

dV4 ■ 0 
dt 

■ 0 
dt 

is ■ ->3 is 
In matrix form, 
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/') 


or U ■ M 

C*1 “ C® ® 0 I 0 j 


or]C - RX 

NOTE: It is recommended that a tree be chosen such that a nonlinear resistance which 

is a link be not associated with a loop (tie set) containing anotirer nonlinear 
resistance. Likewise, a nonlinear resistance which is a tree branch should 
preferably not be associated with a cut set containing anotlier nonlinear resistance. 

Method of Solution 

From the initial state vector X(0) and the curves of the nonlinear elements, tlie 
initial values of the differentiated variables X'(0) are obtained. Tlie time interval of interest 
is now divided into equal increments At, small enough not to impair tlie accuracy of the 
solution, yet large enough not to cause tlie computer to take too long in reaching that 
solution. Increments of the differentiated variables are obtained: 

AX' (At) » A X(0) At 

The values of tlie differentiated variables at time At are: 

^ (At) - X' (0) + A>^' (At) 

Using tlie curves of the nonlinear elements, one can now obtain the state variable values 
at time At, namely, nonlinear resistances exist, 

then X(^*) Otherwise, obtain 

X"(At) X(At) 


V| 

h 

»J 

is 


15 


Utinf tlte curves of the nonlinear resi«tances, a furtlier update of tlie state variables 
% (tit) is obtained, namely, X (At). 

X (At) is now tajken as the initial state vector and tlie same procedure is repeated 
to obtain X (2At), and so on. 

The outputs at any time KAt are given by Y (KAt) ■ B X (icAt) where Y is 
the output vector and B the transmission matrix. Each output is a linear combination 
of the state variables. Tlius, each row of B, say row i, contains the cpefflcients relating 
output i to tlie state variables. 


Computer Program - Flow Chart 
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The program wu implemented on a Hewlett-Packard 9830 computer, provided witlt 
a 9866 printer and a 9862 plotter, A liitini of the program may be found in Appendix A. 

The COMMON itatement permiti modification of individual pieces of data once the 
program terminates. The modified data may be stored on tape witli STORE DATA file# 
command. 

Following is a list of tiie main variables used in the program, togetlier with a description 
of their role. 

C Coefficient matrix for tlie differential equations, 10 x 10. Each row of C multi- 
plied by state-variable vector Z p^uces the derivative of a differentiated variable. 

D Vector of differentiated variables, 10 x I. 

E Output vector at time considered, 4x1. 

F Vector Of non-state variables of dissipative elements, 3x1. 

G Output matrix, 4 x 181. Each row gives the consecutive values of an output in 

time, starting witlr t * 0, 

H Nonlinear characteristics matrix, 10 x 10. Each pair of columns represents a curve; 
hence, up to 5 curves. The first column of a pair contains the x-coordinate values, 
the second column the y-coordinates; for an inductance i and X, for a capacitance 
V and q, for a resistance I and v. 

L Matrix of nonlinear energy-storing element position in Z and corresponding curve 
number in H, 10 x 2. 

N(l) Number of state variables, max. 10. 

N(2) Number of basic time intervals. 

N(3) Print interval. 

N(4) Number of outputs, max. 4. 

N(S) Number of output values to be printed out. 

a- 

N(6) Number of different curves of nonlinear characteristics, 0-S. 
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N(7) Number of nonlliiMr tfMny*itoiing elements, MO. 

N(8) Buie time interval. 

N(9) 

throuf^ Initial condition valuw. 

N(18) 

N(19) Number of nonlinear dissipative elements, 0>3, 

P Matrix for nonlinear dissipative elements, 1x9. To each element corresponds a 

triplet: state variable position in Z, curve number, and kind of state variable (I for 
current, 2 for voltage). 

Q Vector of change in differentiated variables, 10 x 1. 

Q3 Data file number. 

Q4 Number of energy^etoring elements, i.e. number of differential equations. 

R Coefficient matrix for the linear equations of the nonlinear dissipative elements, 

3 X 10. 

T Transmissions matrix. 4 x iO. Each row of T multiplied by state«variabie vector Z 
produces an output value. 

Z State*variable vector at any time, 10 x 1. 

Computer Program - Instructions for Use 

1. NEW JOB (Y OR N)? 

Enter N if you already have a data file for the job; Y if you wish to produce new 
data, and then proceed to question 3. 

2. RLE # ? 

Enter number of existing data file. Proceed to question 20. 

3. NO. OF STATE VARIABLES (1 TO 10)? 

Enter number of state variables. This number may be larger than the order of 
the system, the latter corresponding to the number of fint<order differential 
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aquationf. Ei^ nonlinear diaaipative eJement (resistance, friction) requires a 
state variable, but jives rise to an algebraic equation only. 

4. NO. OF NONLINEAR DISSIPATIVE ELEMENTS (0, I, 2, 3)? 

Enter number of nonlinear dissiipative elements (resistances), up to tliree. 

5. NO. OF NONUNEAR ENERGY-STORING ELEMENTS (0 TO 10)? 

Enter number of nonlinear enerfy-storing elements (inductances, capacitances), 
up to 10. 

6. BASIC TIME INTERVAL (SECONDS)? 

Enter the elementary time interval of integration. For example, if tlie period of 
interest is 1 millisecond and this period is divided into 500 basic time intervals, 
the length of each basic time interval is 2 microseconds, Too large a basic time 
interval speeds up the solution, but impairs its accuracy, 

7. NO, OF TIME INTERVALS? 

Enter number of basic intervals into which tlie total time of interest is divided. 

8. PRINT INTERVAL (I, 2, 3, ...)? 

If tire basic time interval has been chosen judiciously, the variation of the solution 
from one basic time interval to tire next is, in general, too small to be relevant. 
When the solution is eventually printed out, you may want to get tire solution 
every k basic time intervals. This k is called '‘print interval.” For example, let 
the number of basic time intervals be 500. If the print interval is selected as 5 
only solution values # 1, 6, 11, 1^, ... will eventually be printed out, tlrat is 
101 values: the initial condition and 100 intervals. Value # 1 is the initial one, 
before any integration takes place. Incidentally, only those 101 values are kept 
in memory, to be later stored witli the data. Up to 181 output values may thus 
be stored, i.e. NO. OF TIME INTERVALS/PRINT INTERVAL < 180. 
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9. NO. OF OUTPUTS (I, 2, 3, OR 4)? 

&«ttr numbif of output ftinctionii up to 4, An output function ii a linfar combi* 

% 

ration of tli« ittte varUblei. 

SO. TRANSMISSION VECTOR I 

J takM the values of I to the number of outputs desifnated in instruction # 9. 

For output function S, enter tlie coefficients of its linear combination, one at a 
time as requested on display. 

U. COEFFICIENT MATRIX FOR DIFFERENTIAL EQUATIONS 

Enter the coefficients, one at a time, as requested on display. Tliis matrix Is not 
square if noniinear disidpative elements exist. 

’ 12. COEFFICIENT MATRIX FOR ALGEBRAIC EQS. OF N.L. DISSIPATIVE ELEMENTS 
Enter tlte coefficients, one at a time, as requested on display. This information is 
not requested if the answer to instruction 4 is 0. 

13. INIT. COND. VECTOR 

Enter the initial values of all state variables, one at a time, as requested, Tlie 
initial values of the state variables corresponding to nonlinear dissipative elements 
may have to be estimated. Once the computation starts, tlie program will try to 
improve the accuracy of these values to witltin 0.1% in 10 iterations. If it does 
not succeed, it wUl request improved initial conditions and terminate. 

14. NO. OF NONUNEAR CHARACTERISTICS (0 TO 5) ? 

Enter the number of curves describing the characteristics of nonlinear elements. 

This is less than, or equal to, tlie number of nonlinear elements as two or more 
such elements may possess the same characteristics. 

15. ENTER CURVE # J (BY PAIR OF VALUES) 

J takes the value 1 to the number of curves designated in instruction # 14. 

For curve J, enter two values at a time, as requested. The fint value is tlie 
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X'coofdtnitt: Th* currant I for i rcsiitanct ind in inductanct, the voltaie V for ■ 
ctpadtuict. Hie Mcond value ia the y-ooordlnate: the volta«e V for a reiistance, the 
flux linkafe X (volt'oec) for an inductance, and the charie Q for a capacitance. Ten 
Mich pain muit be entered for each curve. This information is not requested if the 
answer to question # 14 is zero. 

16. FOR EACH NONUNEAR ENERGY-STORING ELEMENT, ENTER STATE- 
VARIABLE POSITION AND CURVE # 

For each of the nonUiwar energy-storing elements numbering the answer to question 

# S, enter a pair of values at a time: the state-variable position number in the 
state-variable vector at* entered, for example, under # 13 above, and the associated 
curve number as entered under I S above. More than one element may be associated 
with some curve. 

17. FOR EACH NONLINEAR DISSIPATIVE ELEMENT, ENTER STATE-VARIABLE 
POSITION, CURVE # , AND KIND OF STATE VARIABLE (1 for I, 2 for V) 

For each of the nonlinear dissipative eiements numbering tlie answer to question 

# 4 above, enter a triplet of values at a time: the state variable position number 
in the state variable vector, the associated curve number, and a designation of 1 
if the element is a link or 2 if the element is a tree branch. 

18. STORE DATA (Y OR N) ? 

Enter Y if you wisli all the above data stored on tape, N if you don’t. In the 
latter case skip to instruction # 20. 

19. FILE # ? 

Enter the number of the file where tlie data is to be stored. NOTE: The file 
must have been previously marked (see HP manual tor FIND and MARK commands). 

20. COMPUTE (Y OR N) ? 

Enter V if you wish tlie computation to be carried out in order to obtain solution, 
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N if you don*t. Th« purpoM of Uiii initnictlofi it to Mptnn Hit ditt input from 
tht dita proccttini. Thus, tht wholt Job nt«d not bt dont in ont littini. the 
initnietion it also found useAii in Mpentini the computation job from the output 
Job. 

21. INTERVAL # J. BE PATIENT ! ! ! 

Thii mesMie occun eeery SO intcpation intervals. Nothing to answer. 

22. FINAL STATUS OF STATE-VARIABLE VECTOR IS: 

The solution has ended. The last values of the state variables are printed out. 

They may later be used as initial values to a subsequent time interval. The solution 
and data are stored in the file specified under Instruction # 19 if the answer to 
# 18 is yes. This information may later be retrieved for recomputation, and/or 
printout, and/or plot. Notiting to answer. 

23. TOTAL TIME INTERVAL • - SECONDS 
Notliing to answer 

24. PRINTOUT NEEDED (Y OR N) ? 

Enter Y if you wish to have tire outputs printed out, N if you don’t. 

25. PLOT NEEDED (Y OR N) ? 

Enter Y if you wish to have tin ‘)utputs plotted, and get plotter ready. Enter N 
if you do not need plot, and skip to # 28. 

26. SET PLOTTER, PRESS CONT 

Raise pen, place paper on plotter board, and set all margins. Then press CONT and 
EliECUTE 

27. MIN AND MAX VALUES ? 

Enter pair. It establislies the range of the y*axis. Choose it by referring to the 
values printed out under # 24 above. (MAX>MIN)/S should be a convenient 
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valiM for a marked interval on the y*axis. If the output juat plotted is not the last 
one, the program goes back to*# 26. 

28. END OF JOB 

Nothing to answer. The cassette tape is rewound. 

Possible messages which are followed by program termination; 

1) OUT OF CHARACTERISTICS BOUNDS, PROGRAM ABORTED. 

Variable of nonlinear element is outside range of curve. 

2) PROGRAM ABORTED. IMPROVE INITIAL CONDITIONS. 

, *■ 

Initial values entered for state variables of dissipative elements are not correct. 

I 

3) NO. OF N.L. ELEMENTS - - - NO CURVE. 

Answers to instructions 4 and S are inconsistent with answer to # 14. 

It is recommended tltat tlie data entered be checked for accuracy before computation. 
After entering all data for a new job, answer N to question #20 and all succeeding 
questions. This ends the program. Next, execute tlte following commands to obtain a 
printout of the data entered: 

MAT PRINT N [EXECUTE] 

MAT PRINT T [EXECUTE] 

MAT PRINT C [EXECUTE] 

MAT PRINT R [EXECUTE] 

MAT PRINT H [EXECUTE] 

MAT PRINT L [EXECUTE] 

MAT PRINT P [EXECUTE] 

If corrections are to be made, tlie program does not have to be rerun. Simply enter 
the corrections, one at a time, as follows. If element C (3,5), say, 1$ to be corrected then: 
C (3,5) - new value [EXECUTEJ_ 

and so on. After all corrections are made, store the data. Thus, if tlie data had been 
stored by the program in file # 4, then: 
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STORE DATA 4 


[EXECUTE] 


Now, in order to compute, run tlw profrun again considering tliat you are dealing 
with an old job, 

1. NEW JOB (Y OR N) ? 

N 

2 FILE# ? 

4 

3. COMPUTE (Y OR N) 

V 

Hiis feature is particularly convenient when the circuit model changes due to tlie 
presence of diodes, SCRs, The program would first be run for, say, SO time intervals 
with a diode current as one of tlie outputs. Check the output values printed out. Suppose 
that at interval # 38 the diode current reverses direction. Tlie solution is therefore valid 
only to interval # 37. Then, 

N(2) ■ 37 [EXECUTE] 

N(5) - 38 [EXECUTE] 

STORE DATA same file no. [EXECUTE] 

Rerun the program. At the end, the final values of tlie state variables are to be used 
as initial values for the next run with different equations for a new model and, say, for 
80 time intervals. Then, 

N(2) -80 

N(5) - 81 

N(9) > Z (1,1) and likewise for the rest of the initial conditions. 

Changes in the model 
STORE DATA new file no. 
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Th« solution over the whole time interval of interest may thus be formed of a number 
of solutions over consecutive subintervals. In order to obtairi a single plot for the overall 
solution, tlte next task is to splice sdl subsolutions together, with the proper print 
interval, and store the result in a single file. This can be achieved with the program 
listed in Appendix B and whose details follow. 

1. NO, OF OUTPUTS? 

Enter the number of output functions stored in tlie tiles to be spliced. This number 
must be the same for both files. 

2. TOTAL NO. OF INTERVALS (BOTH FILES)? 

Enter a + b, where a and b are the number of stored data intervals in tiles 1 and 2 , 
respectively. 

3. SUPPLEMENTARY PRINT INTERVAL? 

Both files 1 and 2 to be spliced must also be based on the same print interval, say c. 
If a value d is entered here, it will be considered as a print interval on the output 
data of tlie pre*spUced files. The resulting file will have a print interval of c times d. 

4. FIRST FILE #? 

Enter number of front file. 

5. SECOND FILE #? 

Enter number of rear file. 

6. NO. OF TIME INTERVAIi* LEFT OUT » - - - 

Nothing to answer. It indicates the number of data pieces lost from the end of the 
second file due to the supplementary print interval. For example, let the number 
of intervals of the flnt and second file be 37 and 80, respectively. If the supple* 
mentary print interval chosen is 2, the 80tii output value of tlie second file is lost 
from the spliced combination. 

7. HLE # FOR. NEW DATA ? 

Enter number of file where resulting spliced data are to be stored. 
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10 COM Nn9]»TC4i 10 ]>CC10i 103»ZC 10> 1 If GC4f 1313>HC 10* 10 ]>Ln0i2}*RC3» 101»PC 1 
20 DIM 0C10fl]f£C4fi ]fCtC10il]fFC3fn 
30 PRINT "NEW JOB <V OR N>" 

40 INPUT a* 

90 IF Qt»"r THEN 100 
<0 PRINT "FILE #*• 

70 INPUT 03 
80 LOAD DATA 03 
90 GOTO 1120 

100 PRINT "NO. OF STATE VARIABLES <1 TO lO)" 

110 INPUT Ncn 

120 PRINT "NO. OF NONLINEAR DISSIPATIVE ELEMENTS <0flf2»3>" 

130 INPUT NC193 

140 PRINT "NO. OF NONLINEAR ENERGV-STORING ELEMENTS <0 TO 10)" 

190 INPUT NC73 

180 PRINT “BASIC TIME INTERVAL C SECONDS)" 

170 INPUT NC3] 

130 PRINT "NO. OF TIME INTERVALS" 

190 INPUT NC23 

200 PRINT "PRINT INTERVAL <lf2f3i...)" 

210 INPUT NC33 

220 PRINT "NO. OF OUTPUTS <lf2f3fOR 4)" 

230 INPUT NC43 

240 MAT T«ZERCNC4 3fNm] 

290 Q4>NCn-NC19] 

280 MAT C«ZERC04fNCn3 
270 MAT Z-ZERCNCllf 13 
280 IF NC193-0 THEN 310 
'290 MAT R»ZERCNC193iNn]3 
300 MAT P»ZERClf3#NC1933 
310 NC9 3»INT<NC2 3/NC3 3) + 1 
320 MAT G»ZERCNC4 3tNC533 
330 FOR J-1 TO NC43 
340 PRINT "TRANSMISSION VECTOR" J 
390 FOR I>1 TO NCI 3 
380 DISP "<"n;")"5 
370 INPUT TCJfI3 
380 NEXT I 
390 NEXT J 

400 PRINT- "COEFFICIENT MATRIX FOR DIFFERENTIAL EQUATIONS" 

410 FOR I»1 TO Q4 

420 FOR J»1 TO NCI 3 

430 DISP "<"JIf "f "5J5")"j 

440 INPUT CCIfJ3 

490 NEXT J 

480 NEXT I 

470 IF NCI 9 3-0 THEN 590 

480 PRINT "COEFFICIENT MATRIX FOR ALGEBRAIC EQS. OF N.L. DISSIPATIVE ELEMENTS" 

490 FOR I-l TO NC193 ^ 

900 FOR J«1 TO NCI 3 

9l0 DISP "<"SIJ"f ";j;")"J 

920 INPUT Rn» jn 

930 NEXT J 

940 NEXT I 

990 PRINT "INIT. COND. .VECTOR" 

980 FOR I«1 TO NCI 3 
970 DISP "<"111“)"; 

980 INPUT ZCIfl3 
990 NC3-t>I3-ZCIfl3 
800 NEXT I 


27 



0RK3WAU PA® B 
OF POOR QUALITTf 


6201 

630 

646 

690 

660 

670 

630 

690 

700 

710 

720 

730 

740 

790 

760 

770 

780 

790 

800 

810 

820 

330 

840 

390 

860 

370 

880 

890 

900 

910 

920 

930 

940 

990 

960 

970 

980 

990 

1000 

1010 

1020 

1030 

1040 

1090 

1060 

1070 

1080 

1090 

1100 

1110 

1120 

1130 

1140 

1190 

1160 

1170 

1180 

1190 

1200 


FOR I-NC13+1 TO 10 

Ncs+r>0 • 

NEXT I 

PRINT "NO. OF NONLINEPR CHftRflCTERISTICS <0 T0-9>" 

INPUT NC6] 

IF NC 61*0 THEN 970 
REDIM Hn0>2#NC6]] 

FOR J*lTO NC63 

PRINT "ENTER CURVE #"J5"<BY PRIR OF VfiLUES)" 

FOR 1*1 TO 10 
niSP "PAIR #"?I! 

INPUT HCI»2*J-1]>HCIj2*J] 

NEXT I 
NEXT J 

IF NC7 3-0 THEN 890 
REDIM LCHC73>2] 

PRINT 

PRINT "FOR ERCH NONLINERR ENERGY-STORING ELEMENT?" 

PRINT "ENTER STRTE-VflRIRBLE POSITION AND CURVE #" 

FOR 1*1 TO NC71 
DISP "<"?i;”>"? 

INPUT UCI?l]?LCI>2] 

NEXT I 

GOTO 860 

MAT L»ZERC1? 1 3 

IF NCI 9 3*0 THEN 940 

PRINT 

PRINT "FOR EACH NONLINEAR DISSIPATIVE ELEMENT? ENTER STATE-VARIABLE" 
PRINT "POSITION? CURVE #? AND KIND OF STATE VARIABLE <1 FOR I? 2 FOR V)' 
FOR I»1 TO NU9 3 

INPUT PC 1?3#I-23?PC lf3irl-l 3?PC 1?3»1 3 

NEXT I 

GOTO 1050 

MAT P-ZERC I ? 1 3 

MAT R*ZERC1?13 

GOTO 1090 

IF NC7 3+NC19 3-0 THEN 1O10 

PRINT "NO. OF N.L. ELEMENTS*" ? NC 7 3+NC 19 3? " ? NO CURVE." 

DISP "ERROR! PROGRAM ABORTED." 

STOP 

MAT P-ZERC 1? 1 3 
MAT R-ZERC I? 1 3 
MAT H=ZERC1?13 
MAT L»ZERC1?13 

PRINT "STORE DATA a OR N>"? 

INPUT Rf 

IF RI»"N" THEN 1270 
PRINT "FILE #"? 

INPUT 03 
STORE DATA 03 
GOTO 1270 
04*NC 1 3-NC 193 
R*«"Y" 

REDIMTnrC4 3?NC 1 3 3? CC04?NC 1 3 3?ZCNC 1 3? 1 3? GC NC 4 3? NC 9 3 3 

IF NC63-0 THEN 1260 

REDIM HC10?2^NC6 33 

IF NC73»0 THEN 1200 

REDIM LCNC73?23 

GOTO 1210 

REDIM LC1?13 
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1210 IF NC19I-0 THEN 1240 

1220 REDIH RCNM93»NCn]»PCii3»NCl9]3 

1230 GOTO 1270 

1240 REDIH RClil3rPCl»13»FCM3 
1250 GOTO 1270 

12€0 REDIH HCt»13>Ln!l3»RClrl3rPChl3 
1 270 RED IH DC 04 » 1 3f EC NC 4 3 f 1 3 > QC 04 f 1 3 
1200 IF NU9 3-0 THEN 1300 
1200 REDIH FCHC10 3>13 
1300 PRINT “COMPUTE <Y OR N)“ 

1310 INPUT Of 

1320 IF Ot*"N“ THEN 2040 

1330 FOR I»1 TO NCI 3 

1340 sci^naNco-t-n 

1350 IF I >04 THEN 1370 

1360 DC1»13»NC8+I3 ^ 

1370 next I 

1330 IF NC63-0 THEN 1640 
1390 IF NC73-0 THEN 1460 
1400 FOR Jl-1 TO NC73 
1410 Pl«DCLCJl»13j 1 3 
1420 P3-LCJlf23 
1430 GOSUB 2730 
1440 DCLCJ1 j13»13»P2 
1450 NEXT J1 

1460 IF NC193«0 THEN 1640 
1470 FOR Qlsi TO 10 
1480 HfiT F-R*2 
1490 02-0 

1500 FOR Jl-1 TO NC19*3 

1510 Pl-FCJl* 1 3 

1520 P3»PCl»3i^Jl-l 3 

1530 GOSUB PCl»3*J13 OF 2900 *2730 

1540 IF 2CQ4+Jlrl3»0 THEN 1560 

1550 02»Q2+flBS<<2CQ4+Jlj 13-P2>/ZCG!4+Jl»i 3) 

1560 2C94'*‘Jlil 3»P2 
1570 NEXT J1 

1530 IF O2<0.OO1»NC193 THEN 1620 
1590 IF 01410 THEN 1630 

1600 PRINT "PROGRRM RBORTED. IMPROVE INITIRL CONDITIONS.” 
1610 STOP 
1620 Ql»l0 
1630 NEXT 01 

1640 MRT E»T*2 , 

1650 FOR I»1 TO NC43 
1660 GCMl-ECIfl] 

1670 NEXT I 

1630 FOR Kl-1 TO NC23 

1690 MfiT 0«C*Z 

1700 MfiT 0-<NC8 3)#0 

1710 MfiT D-D+0 

1720 REDIM ZC04rl3 

1730 MfiT Z-0 r 

1740 REDIM ZCNC13il3 

1750 If NC63-0' THEN 1950 

1760 IF NC 7 3-0' THEN 1S30 

1770 FOR Jl-1 TO NC73 

1730 Pl»ZCUJl>13rl3 

1790 P3-LCJ1>23 

1300 GOSUB 2900 
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1810 zcLCJi>nai-P2 
1820 NEXT Jl 

1830 IF NCr?>0 THEM 1830 

1340 FOR Ol-l TO 2 

1850 MfiT F-R»2 

i860 FOR Jl«l TO HC193 

1870 Pl»FUl»n 

1880 P3«PC1>3»J1-1I 

1890 GOSUB PClf3»J13 OF 2900f273O 

1900 ZCQ4+Jl>n»P2 

1910 NEXT J1 

1920 NEXT Q1 

1930 IF K14aHT<Kl/50>>*50 THEN 1950 
1940 PRINT "INTERVfiL 4“Kr*BE PRTIENT.!!!" 

1950 IF K1#<INT<K1/NC33>>#NC3] THEN 201O 
1960 MfiT E»T*2 
1970 K2»Kl/NC3]+l 
1980 FOR I»1 TO NC4] 

1990 GCIrK23«EC If13 
2000 NEXT I 
2010 NEXT K1 

2020 IF R#«"N" THEN 2040 
2030 STORE DflTfl Q3 
2040 PRINT 

£050 PRINT "FINAL STATUS OF STATE-VARIABLE VECTOR IS!" 
2060 REDIN ZClrNC in 
2070 MAT PRINT Z 

£080 PRINT "TOTAL TINE 1NTERVAL»"NC2 3*NC 8 3"SEC0NDS" 
£090 PRINT 

2100 DISP "PRINTOUT NEEDED <V OR N)"? 

2110 INPUT 13* 

£120 IF 0*»"N" THEN £210 
2130 PRINT 

2140 FOR I»1 TO NC43 
2150 PRINT "OUTPUT"! 

£160 FOR J«1 TO NC5] 

£170 PRINT GCIjJ].! 

2180 NEXT J 
2190 PRINT 
2200 NEXT I 

2210 DISP "PLOT NEEDED fV OR N.V5 

2220 INPUT Qt 

2230 IF THEN £670 

2240 FOR I-l TO NC43 

£250 DISP "SET PLOTTER? PRESS CONT" 

2260 STOP 
£270 PEN 

£200 DISP "MIN AND MAX VALUES"? 

2290 INPUT E1?E2 

2300 SCALE 0?l.UNC2 3f^NC8 3?El?E2 
?1310 E3»INT<NC2 3/100>#10#NC8 3 
£320 E4-AB84Ea-El) 

2330 E5-LGT(E4> 

£340. E4-E4/10tE5 

2350 ES-rNT<E4)*10t<E5-l) 

2360 IF E2»E1>0 THEN 2390 

2370 E4-0 

2380 GOTO 2430 

2390 IP E1<0 THEN 2420 

2400 E4»E1 
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24 Id, GOTO 2430 
24201 E4-E2 

2430 E5-1.1«NC2]#HC8] 

2440 XflXIS E4iE3»0>E5 v? 

2450 VAXI3 0>E€»Ei?E2 

24E0 LABEL <#> 1.5r 1.7>0»3.5/'11> 

2470 FOR V«0 TO lO STEP 2 

2400 IF A63<El>V»E6-E4KEe THEN 2520 

2490 PLOT 0>Y#E6^Ebl 

2500 CPLOT 2»-0.3 

2510 LABEL <260OXY#E6-»-En 

2520 NEXT Y 

2530 LABEL 1 ,5> 1 . 7»Pl^’2»8.5/l 1 > 

2540 FOR X»0 TO lO STEP 2 
2550 IF X»0 THEN 2590 
2560 PLOT X»E3>0»1 
2570 CPLOT 2>-0.3 
2580 LABEL (260OKX»E3> 

2590 NEXT X 

2600 FORMAT E8. 1 

2610 PLOT 0»GCI»n»-2 

2620 FOR J«2 TO NC53 

2630 PLOT <J-l)#NC31*NC8J5GCI»JJr2 

2640 NEXT -J 

2650 PEN 

2660 NEXT I 

2670 REWIND 

2630 PRINT 

2690 PRINT "END OF JOB" 

2700 END 

2710 REM 8 FROM STATE VAR. TO PROBLEM VAR. 

2720 STOP 

2730 K3«2#P3-I 

2740 K4»2*P3 

2750 FOR I«1 TO 10 

2760 IF Pl»HCIfK3] THEN 2S20 

2770 IF Pl>HnfK33 THEM 2840 

2780 IF I»1 THEN 2850 

2790 J»I-1 

2800 P2«HC J > K4 J+ < HC I K4 3-HC J > K4 ] > 4 < P I -HC J » K3 3 ) / C HC I > K3 3-HC J > K3 3 ) 

2810 GOTO 2370 

2320 P2-HCI»K43 

2330 GOTO 2370 

2840 NEXT I 

2850 PRINT "OUT OF CHARACTERISTICS BOUNDS. PROGRAM ABORTED." 

2360 STOP 
2370 RETURN 

2880 REM 8 FROM PROBLEM VAR. TO STATE VAR. 

2390 STOP 

2900 K3-24P3-1 

2910 K4»2#P3 

2920 FOR I«I TO 10 

2930 IF P1«HCIjK43 THEN 2990 

2940 IF P1>HCI»K43 THEN 3010 

2950 IF 1^1 THEN 3020 

2960 J»I-1 

2970 P2»HCJ»K33+CHn»K'33-H[ J»K33)#CP1-HC JjK 43)/<HC IfK4 3-HC JjK43) 

2980 GOTO 3040 

2990 P2»HCI»K33 

3000 GOTO 3040 

3010 NEXT I 

3020 PRINT "OUT OF CHARACTERISTICS BOUNDS. PROGRAM ABORTED." 

3030 STOP 
3040 RETURN 
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to COH Hn91>TC4>L0i>Cn0<101>SCt0>n>CC4itSl]»HCl0»t0MCt0>S]»Rt3tt0nFni91 
20 DIM AC4»101] 

30 DISP "HO. OF OUTPUTS"* 

4i; INPUT E4 

50 PRINT "TOTfiiL HO. OF IHTERVftLS (BOTH FILES)" 

SO INPUT E9 

70 PRINT "SUPPLEHENTPRY PRINT INTERVAL" 

80 INPUT ES 
SO ES*INT(ES/ES)+l * 
too REDIM RCE4 fE3] ‘ 

110 DISP "FIRST FILE 
120 INPUT 03 
130 LOAD DATA 03 
t40 IF E4-NI4] THEN ISO 

150 PRINT "ERROR! DIFFERENT HO. OF OUTPUTS." 
tS0 DISP "PROGRAM TERMINATED" 

170 STOP 

130 E3»NI33 

130 E2-NC23 

200 E5-NI53 

2t0 REDIM CC£4iE51 

220 ES«1 

230 FOR J»1 TO £5 STEP E6 

240 FOR I"l TO £-4 

250 ACI*E3 3-GC IrJI 

2S0 NEMT I 

270 £3«E3+l 

230 NEXT J 

230 £7«E5-<(E8-2:>#Eo+l) 

380 E73ES-E7+1 

310 REDIM GC4> 131 1 

320 DISP "SECOND FILE #"T 

330 INPUT 03 

340 LOAD DATA 03 

350 IF E4-NC4I THEM 370 

3«0 GOTO 150 

370 IF E3-HC31 THEN 400 

330 PRINT "ERROR! FILES DIFFER IN PRINT INTERVAL." 

330 GOTO ISO 

400 IF E54NC51-2»E3 THEN 430 

410 PRINT "ERROR! TOTAL NO. OF INTERVALS -"(ES+NCSI-S) 

420 GOTO ISO 

430 REDIM GCE4,MC5n 

440 FOR J-E7 TO NC51 STEP E3 

450 FOR I»l TO E4 

460 AC I>E3>GC I*J] 

470 NEXT I 
430 E3*E3tl 
430 NEXT J 
500 E3-E8>t 

510 E7»(E5+NC5I-2>-<E3-l)*E6 

520 PRINT "NU. OF TIME INTERVALS LEFT 0UT»"JE7 

530 REDIM GCE4.ES1 

540 MAT G«A 

550 MC51-E8 

560 HC2I-NI2HE2 

570 HC33»NC3 3*eS 

530 DISP "FILE i FOR HEW DATA"? 

530 INPUT 03 

600 DISP "IF TAPE READYf PRESS CONT" 

610 STOP 

620 STORE DATA 03 

630 DISP "PRINTOUT NEEDED (Y OR N)"5 

640 INPUT 0» 

650 IF 0*»"N" THEN 670 

660 MAT PRINT G 

670 DISP "END OF JOB" 

630 END 
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